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In this work we investigate particle acceleration in supernova remnant shocks within a semi- 
analytical formalism which self-consistently accounts for particle acceleration, amplification of 
the magnetic field via streaming instability and back-reaction of both accelerated particles and 
magnetic turbulence on the shock dynamics. In particular, we study the interplay between parti- 
cle injection and magnetic field amplification, showing how a phenomenological but reasonable 
saturation of the standard streaming instability may lead to quite steep spectra for the accelerated 
particles. We comment on the implications that such a scenario may have on the comprehension 
of the diffuse spectrum of Galactic cosmic rays and of gamma-ray observations of single SNRs. 
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1. The standard scenario 

The association between supernova remnants (SNRs) and Galactic cosmic rays (CRs) has been 
very popular since 1934, when Baade and Zwicky argued that this class of astrophysical objects can 
account for the required CR energetics [|T]|. Along with this argument, the so-called SNR paradigm 
for the origin of CRs has been supported also by E. Fermi's proposal of a very general mechanism 
for particle acceleration, which turned out to be particularly efficient for particles diffusing around 
SNR forward shocks [Q, |3|]. This diffusive shock acceleration (DSA) mechanism is also very ap- 
pealing because it naturally returns power-law spectra 0(£') °« E^^ for the energy distribution of 
accelerated particles, with F = 2 for large sonic Mach numbers. Moreover, such a spectral index 
does not depend on the details of the scattering processes but only on the compression ratio felt by 
diffusing particles. 

1.1 The diffuse spectrum of Galactic CRs 

By measuring the ratios between secondary and primary CRs, it is possible to infer that CRs 
diffuse during their propagation in the Milky Way and that the Galactic residence time scales as 
E^^, with 5 ~ 0.3 — 0.6. Since the CR spectrum measured at Earth is oc £^2.7 fj-Qjjj about a GeV 
up to the knee (around 3 x lO^GeV), in this wide range of energies we have a strong observational 
constraint, namely F + 5 ~ 2.7 and, therefore, F = 2.1 — 2.4. CR spectra at the sources may thus 
be steeper than the ones predicted by standard DSA at SNR forward shocks 

This point is made even more problematic by the fact that, in order to account for correct ener- 
getics, SNR shocks have to channel large fractions (from 10% to 50%, depending on the Galactic 
propagation details) of their kinetic energy into accelerated particles. Therefore, the SNR paradigm 
for the origin of CRs naturally suggests SNR shocks to be modified by the dynamical back-reaction 
of particles they accelerate: any reliable DSA theory has to be non-linear (NLDSA), i.e. has to in- 
clude the feedback of CR pressure and energy in the shock hydrodynamics. 

The most important output of the kinetic NLDSA theory, independently of the approach 
adopted for describing the coupling between CR transport and SNR hydrodynamics [^], is that the 
spectrum of accelerated particles has to be concave, i.e. steeper (harder) than E^^ — the standard 
prediction of DSA — at low (high) energies, as a consequence of the shock modification induced 
by an efficient particle acceleration ^ |8]|. This prediction leads to some very natural questions: 
why such a concavity is not observed in the diffuse spectrum of Galactic CRs? And also, how can 
NLDSA be consistent with source spectra steeper than E^'^ up to the knee? 

Before trying to answer these questions it is necessary to stress that, in principle, there is 
no one-to-one correlation between the instantaneous spectrum of particles being accelerated at a 
shock, whose slope is predicted by Fermi's mechanism, and the spectrum of particles released 
from a SNR during its lifetime. The latter is in fact the convolution over different evolutionary 
stages of (concave?) spectra, with different normalizations and cut-offs. In particular, during the 
Sedov-Taylor stage the instantaneous maximum energy decreases with time as a consequence of 
the decrease of both shock velocity and magnetization level: at any given time, the highest-energy 
particles cannot but escape the system from upstream |^. A very simple argument demon- 
strates how the self-similar adiabatic evolution of a SNR during its Sedov-Taylor stage may also 
lead to a ~ E^-^ spectrum for these escaping particles. The actual details of how and when acceler- 
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ated particles are released into the Galaxy and become CRs are still poorly understood. Different 



possibilities for particle escape are investigated in Ref. []10|], where it is also shown that, most likely, 
escape from upstream is necessary to account for the highest -energy part of the broadband spec- 
trum, while at lower energies, say below 1 TeV, the source spectrum is almost invariably dominated 
by the accelerated particles which have been advected downstream, suffered adiabatic losses, and 
finally released at the "death" of the remnant itself. 



In all of the NLDSA calculations of the most abundant CR species [11, 12, 13] the spectral 
slope at the source has never been found to be larger than r~2.1— 2.15,in turn implying 5 > 0.55. 
Nevertheless, as discussed e.g. in §7 of Ref. [Q], a strong energy dependence in the residence time 
is at odds with the low level of anisotropy measured in the CR arrival direction, which seems to 
point towards the lowest limit consistent with the measured B/C ratio, namely d ~ 0.3. 

1.2 Hadronic emission from SNRs 

There is another independent piece of information coming from the present generation of y-ray 
telescopes, which are providing us with an unprecedented wealth of data from sources coincident 
with SNRs. Very briefly, there are two main scenarios in principle able to account for the high- 
energy photons detected in the GeV band by Fermi and AGILE and in the TeV band by MAGIC, 
VERITAS, HESS, CANGAROO, MILAGRO, etc.: the leptonic scenario, in which y-rays are due 
to inverse Compton scattering or bremsstrahlung by relativistic electrons, or the hadronic scenario, 
in which they are due to the decay of neutral pions produced in nuclear collisions between the 



relativistic nuclei and the background plasma [14|. Even if there are no clear-cut evidences of the 



predominance of a scenario over the other yet, we cannot but notice that basically all of the SNRs 
detected in /-rays show a photon spectrum softer than E^'^ [p3|]. 

This fact is particularly intriguing because inverse Compton scattering returns rather hard pho- 
ton spectra, while both relativistic bremsstrahlung and pion decay provide a one-to-one map of the 
parent electron and proton spectra in terms of spectral slope: it is very likely that by observing 
7-ray bright SNRs we can probe the instantaneous spectrum of accelerated particles and, in turn, 
have an insight into how NLDSA effectively works. 

In this respect, it is interesting to notice that most of the detected SNRs show spectra in the 
range E^^-^ — E^^-^, in nice agreement with the adoption of a mildly energy-dependent Galactic 
residence time and with the measured anisotropy. Putting all these evidences together, one may 
conclude that the most natural acceleration mechanism at work in SNRs has to produce spectra 
0(£') oc E^^-^, hence significantly steeper than what predicted by standard NLDSA calculations. 



2. CRs and magnetic field amplification 

We want to exploit the semi-analytical formalism put forward in Refs. [16, [T^ to investigate 
the possibility to obtain such a steep source spectrum in the context of a NLDSA theory by in- 
cluding also the effects induced by a very efficient magnetic field amplification (MFA). The strong 
non- thermal X-ray luminosity of young SNRs has been in fact explained as due to synchrotron 
emission by relativistic electrons radiating in magnetic fields as high as a few hundred /iG, i.e. 
orders of magnitudes larger than the typical interstellar ones (e.g. [17] and references therein). 
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These strong magnetic fields are very likely produced by the same accelerated particles via 
plasma instabilities which, in the quasi-linear regime, saturate when the ratio between the pressure 
in accelerated particles and in magnetic turbulence is of order Pct/Pb — 2Ma, where Ma = VsIi/Va 
is the Alfvenic Mach number of the shock We extrapolate this result to the non-linear regime 
5B » Bo by assuming that the relevant Alfven velocity has to be calculated in the amplified mag- 
netic field 5B, and not in the background one as prescribed by the quasi-linear theory. Additional 



comments and caveats related to this choice can be found e.g. in Refs. [ ]15[ , |18| ]. 

Very generally, there are several effects induced by a very efficient MFA: 1) the pressure in 
magnetic turbulence may easily become dominant with respect to the thermal plasma pressure. 



governing the hydrodynamics and preventing an excessive shock modification |[18|, [19|]; 2) the 
velocity of the waves particles scatter against, which is of order ~ Va, may become a non-negligible 
fraction of the fluid velocity, altering in a significant way the compression ratios actually felt by CRs 
and eventually making their spectrum steeper than E^^ [Q, |l^, 15]; 3) since Ma{5B) <^ Ma (Bo), 
for a given gradient in the CR pressure the instability saturation might be achieved for larger SB. 

We consider a SN explosion releasing 1.4 solar masses and 10^1 erg in the interstellar medium 
with density, temperature and magnetic field equal to 0.01 protons cm^^, lO^K and 5/iG, as in 



Ref. [10]. The self-similar shock evolution is coupled with a non-linear treatment of particle ac- 
celeration, shock dynamics and MFA, including particle injection from the thermal bath [pO|]. In 
particular, all the protons with momentum i^inj times larger than the thermal momentum down- 
stream are assumed to be injected into the acceleration process (thermal leakage): the smaller c^inj, 
the larger the fraction of injected particles. 

MFA amplification upstream is taken as due to resonant streaming instability, producing a 
magnetic field given by 

5B^ Per 

=Pb^^^. (2.1) 

871 2Ma 

The dynamical back-reaction of such an amplified magnetic field on the shock dynamics is also 



included [lO, 18]. Particle escape due to the decrease of the SNR confining power during the 



Sedov-Taylor stage is accounted for by imposing the CR distribution function to vanish beyond a 



free-escape boundary placed at a distance O.SRsh upstream of the shock |16||, while escape from 
downstream is treated as in Ref. [ |T0| ] by posing j3 =0.1. 

In the left panels of Fig. |l] the total CR spectra produced by our benchmark remnant are shown 
(thick solid lines) as the sum of three different contributions: escape from upstream (dashed lines), 
leakage from downstream (dotted) and release at the end of the Sedov-Taylor stage, when the SNR 
non-thermal activity is expected to fade (dot-dashed). In the right panels, instead, some snapshots 
of instantaneous proton spectra are shown; the first age corresponds to the transition from the 
ejecta-dominated to the adiabatic stage, when energies as high as the knee are achieved [p|]. 

We consider here two different configurations for the velocity of the scattering centers, in order 
to illustrate how crucial this ingredient may be. The top panels of Fig. |I|show a case in which MFA 
is at work but scattering centers comove with the fluid, i.e. Va = 0. This case is compared with the 
panels below, where Va is included and calculated in the amplified magnetic field 5B, instead. 

There are two very important points to be noticed. First, considering an effective Alfvenic 
drift in the amplified magnetic field leads to spectra much steeper than in the case Va = 0, where 
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Figure 1: Source spectra without (Va — 0, top panels) and with the Alfvenic drift included (Va{8B), other 
panels), for different injection efficiencies — 3.9, 3.6, 3.3, as in the legend. The solid lines in the left pan- 
els represent the total CR spectrum, distinguished into contributions due to escape from upstream (dashed), 
leakage from downstream (dotted) and release at the SNR death (dot-dashed lines). In the right panels 
instantaneous spectra of accelerated panicles are instead shown as a function of the SNR age. 
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the convolution over time of concave spectra leads to a power-law spectrum slightly flatter than 
E^^, as already shown in Ref. [p^]. 

Second, when the Alfven velocity (and hence the saturation of the resonant streaming instabil- 
ity) are taken in the amplified magnetic field, total spectra are straight power-laws whose spectral 
indexes correlate with the parameter which tunes the thermal leakage efficiency, i^inj , here taken to 
vary between 3.3 and 3.9. In particular, second, third and fourth rows in Fig. |]have increasingly 
larger thermal leakages, corresponding to inject into the acceleration process fractions as large 
as ~ 10"^ ((^inj = 3.9), ~ 10"^ ((^inj = 3.6) and ~ lO"'^ (i^inj = 3.3) of the particles crossing the 
shock. We can interpret this behavior as follows: the larger the number of particles injected, the 
larger the CR pressure and, as a consequence of Eq. 2A , the larger the self-generated magnetic 
turbulence. This last process, however, produces a crucial non-linear feedback which makes the 
particle spectrum steeper and eventually reduces the amount of energy channelled into CRs. The 
balance between these two opposite effects returns self-consistent solutions with both efficient par- 
ticle acceleration (around 10%) and, occasionly, very steep spectra, like the case i^inj = 3.3 in Fig. |l|. 
The spectral slopes of the total SNR spectra are in fact F ~ 2.15,2.28 and 2.45 for the cases with 
Va{SB). The fine structure of such a non-linear interplay will be widely discussed in a forthcoming 
paper. 



3. Conclusions 

The details of how MFA occurs in SNR shocks may radically change our previous comprehen- 
sion of efficient DSA theories, leading to a revised interpretation of many observational facts. In 
order to illustrate this point, we calculated instantaneous and total CR spectra produced at a typical 
SNR forward shock by adopting a rigorous semi-analytical kinetic approach to NLDSA, on top of 
which we added a phenomenological description of very efficient Alfvenic turbulence generation 
due to resonant streaming instability. 

We showed that it is possible to account both for quite large acceleration efficiencies (around 
10%) and, at the same time, for CR spectra significantly steeper than what predicted by previous 
NLDSA calculations. More precisely, when the Alfven velocity, Va, is calculated in the amplified 
magnetic field, the intrinsic non-linear interplay between accelerated particles, magnetic field am- 
plification and shock dynamics provides steeper and steeper spectra for larger and larger number 
of injected particles, actually reversing the trend one would obtain by taking Va in the background 
magnetic field, as prescribed by quasi-linear theory. 

In terms of slope, normalization and cut-off the spectra we obtained are consitent with basi- 
cally all the observational evidences inferred from the diffuse spectrum of Galactic CRs and from 
7-ray-brigth SNRs [||]. 
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